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In tunneling spectroscopy studies of ferromagnet/superconductor (F/S) junctions, the effects of 

O . spin polarization, Fermi wavevector mismatch (FWM) between the F and S regions, and interfacial 

^ ' resistance play a crucial role. We study the low bias conductance spectrum of these junctions, gov- 

pL^ , erned by Andreev reflection at the F/S interface. We consider both d- and s-wave superconductors 

as well as mixed states of the d + is form. We present results for a range of values of the relevant 

parameters and find that a rich variety of features appears, depending on pairing state and other 

conditions. We show that in the presence of FWM, spin polarization can enhance Andreev reflection 

C ' and give rise to a zero bias conductance peak for an s-wave superconductor. 

y ; 74.80.Fp, 74.50+r, 74.72-h 
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^ ■ I. INTRODUCTION 

-)— > ■ 

<^ 

The development and refinement in recent years of new techniques in materials growth has made it possible to 

I , fabricate superconducting heterostructures with various materials and high quality interfaces. Thxisa advances, coupled 

■ O ' with the continuiag intense level of activity in the study of the nature of high temperatura^^B and other exotic 

^ \ superconductors,QLl has led ta.fenewed interest in tunneling spectroscopy. 

jx ■ It has been demonstratedaQD that this technique yields information about both the magnitude and the phase of 

1—^' the superconducting pair potential (PP). This implies that the method can provide a systematic Sf^y to distinguish 

among various proposed PP candidates, including both spin si|H||let pad spin triplet pairing states.™ For example, it 

V^ ■ has been argued that the observed zero bias conductance peakQu'EErO (ZBCP), attributed to mid-gap surface states, 

is an indication of unconventional superconductivity with a sign change of the PP, as it occurs in pairing with a 



> 

QQ ' dx2_y2-waMe symmetry. Furthermore, the splitting of the ZBCE-and the forming of a finite bias peak (FBCP) in 



f~^ I the conductance spectrum has been examined and interpretedcJ iB as support for the admipslurfi of an imaginary PP 



CNJ ■ component to the dominant dx2-y2-wa.ve part, leading to a broken time-reversal symmetryJljO 

^D I The same developments, and the ability to make low interfacial resistance junctions between high spin polarization 

O^ , ferromagnets and superconductors, have stimulated siffiiificant efforts to study transport in these structures .Il3 There 

^/ ' have been various experiments in both-jCapventionaLjiill and high temperature superconductor£3~Ej (HTSC's), as 

well as re-examinations of earlier workn3'EZl which was performed generally in the tunneling limit of strong interfacial 

barrier. Theoretical studies of the effects of spin polarized transport on the current-voltage charActeristics and the 

conductance in ferromagnet/supercanductor (F/S) junctions have been carried out in conventionaEs and, recently, in 

high-temperature superconductors.EEnEl The feasibility of nanofabricating F/S-|Structures has also generated interest 

f~| ' in studying the influence of ferromagnetism on mesoscopic superconductivityH3 

O ■ One of the important questions raised by the p ossib ility of making high transmissivity F/S junctions was that of 

^ ■ studying the influence of Andreev reflection (AR)c§E3~Ej on spin polarized transport. In AR an electron, belonging 

J> , to one of the two spin bands, incoming from the ferromagnetic region to the F/S interface will be reflected as a hole 

! '~j ' in the opposite spin band. The splitting of spin bands by the exchange energy in ferromagnetic materials implies that 

rS ' only a fraction of the incoming iQCJdent majority spin electrons can be Andreev reflected^ This simple argument 

j^ ■ was used in previous studiesLJcj^ES to infer that the effect of spin polarization (exchange-energy) was generally to 

■ ■ ' reduce AR. The sensitivity of AR to the exchange energy in a ferromagnet was employedlljC3 to determine the degree 

of polarization in various materials. 

In this paper we will study the tunneling spectroscopy of F/S junctions. We will adopt the basic approach of Ref. ^ 
but we will extend and generalize it to include the effects of spin polarization, . the .presence of an unconventional PP 
state (pure or mixed) , and the existence of Fermi wavevector mismatch (FWM)LZtt3 stemming from the different band 
widths in the two junction materials. Our aim in this paper is twofold: Firstly, to investigate and reveal novel features 
in the conductance spectra arising from the interplay of ferromagnetism and unconventional superconductivity, and 
secondly, to show the importance of FWM, and how its inclusion can lead to some unexpected results, even for F/s- 



wave superconductor junctions, where we find, for example, that in some cases Andreev reflection can be enhanced 
by spin polarization. 

In the next section (Sec. O), we present the methods we use to obtain the amplitudes for the various scattering 
processes that occur in the junction when spin polarized electrons are injected from the F into the S region. We will use 
these methods to calculate the conductance of the F/S junctions. In Sec. |lD, we first give results for a conventional (s- 
wave) superconductor in the S side, and then illustrate the unconventional case of the pairing potential by considering 
both pure d- and mixed d + is-wave symmetry. In Sec. [V, we summarize our results and discuss future problems. 



II. METHODS 



As explained in the Introduction, we investigate in this work F/S junctions by extending and generalizing the 
techniques previously employed in the study of simpler cases without spkii-BaLarpation, or for conventional supercon- 
ductors. Thus, we use here the Bogoliubov-de Gennes (BdG) equationsBuE3'E3'L3 in the ballistic limit. We consider 
a geometry whece the ferromagnetic material is at x < 0, and is described by the Stoner model. We take the 
usual approached of assuming a single particle Hamiltonian with the exchange energy being therefore of the form 
h{r) = hoQ{—x), where Q{x) is a steuiunctipn. The F/S interface is at x = 0, where there is interfacial scattering 
modeled by a potential V{r) — H6{x)ptBlEjE.-3fid H is the variable strength of the potential barrier. The dimension- 
less parameter characterizing barrier strengthtSI is Zq = mH/hkF, where thft-jaffectiyiajniass, m, is taken to be equal 
in the F and S regions. In the superconducting region, at x > 0, we assumeOQ SEj'lJ that there is a pair potential 
A(k',r) = A(k')8(x). This approximation for the PP becomes more accuratec^ in the presence of FWM and allows 
analytic solution of the BdG equations. We will denote quantities pertaining to the S region by primed letters. 

P&om these considerations, the BdG equations for F/S junction, in the absence spin-flip scattering, can be written 
a^ 



Ho - psh A 

A* -{Ho + psh) 



us 



us 



(2.1) 



where Hq is the single particle Hamiltonian and pg = ±1 for spin S =T,i- 
PP A are as defined above. The excitation energy is denoted by e, and us-. 



The exchange energy h{r) and the 
v-^ are the electronlike quasiparticle 



(ELQ) and holelike quasiparticle (HLQ) amplitudes, respectively. We take Ho = — 7i^V^/2m -I- ^(r) — E/ , where 



V{r) is defined above 
Ef = {h^k%J2m 



Ef 



h^e 



E' 



h^k'^ 



In the 

Fi/2m)/2 



F region, we have Ep 



We 



Ep = h kp/2m, so that Ep is the spin averaged value, 
assume that in general it differs from the value in the superconductor. 



p/2m. Thus, we take the Fermi energies to be different in the F and S regions that is, we allow for 



different band widths, stemming from the different carrier densities in the two regions. Indeed, as the results in the 
next Section will show, the Fermi wavevector mismatch (FWM) between the two regions has an important influence 
on our findings. We will parameterize the FWM by the value of Lq, Lo = k'p/kp and describe the degree of spin 
polarization, related to the exchange energy, by the dimcnsionless parameter X = ho/ Ep. _. 

The invariancc of the Hamiltonian with respect to translations parallel to a; = implies conservatiorcS of the (spin 
dependent) parallel component of the the wavevector at the junction. As we shall show, this will be an important 
consideration in understanding the possible scattering processes. An electron injected from the F side, with spin 
S =t,i, excitation energy e, and wavevector k^ (with magnitude ^ipn {2m,/h?Y/'^\Ep -I- e + pshoY^"^), at an 
angle 6 from the interface normal, can undergo four scattering processesQ'tj each described by a different amplitude. 
Assuming specular reffection at the interface, these can be characterized as follows: 1) Andreev reffection, with 
amplitude that we denote by ag, as a hole with spin, 5', belonging to the spin band opposite to that of the incident 
electron {pg = —ps), wavevector k-. {k%- — {2m/h'^y^'^[Ep — e + pgho]^^'^), and spin dependent angle of reflection 

9-g, generally different from 6'.t3 As is the case with the angles corresponding to the other scattering processes, 6^, as 
we shall see below, is determined from the requirement that the parallel component of the wavevector is conserved. 
Even in the abseiice of exchange energy (ho = 0), one has that, for e j^ 0, 9-^ (although then spin independent) is 
slightly differento from 9. When ho > 0, the typical situation is, as we discuss later, that \9j\ < \9\ < \9j\. 2) The 
second process is ordinary reffection into the F region, characterized by an amplitude which we call 65, as an electron 
with variables S, — kj, —9. The other two processes are: 3) Transmission into the S region, with amplitude C5, as 
an ELQ with k'<y, and 4) Transmission as a HLQ with amplitude ds and wavevector — k|^. Here the corresponding 
wavevector magnitudes are k'^ — {2m/h'^Y/'^[E'p ± (e^ — |As±p)^/^]^/^. We denote by As± = |As±| exp(i(/)s±), the 
different PP's felt by the ELQ and the HLQ, respectively, as determined by k^ . We see, therefore, that up to four 
different energy scales of the PR afa involved for each incident angle 9. In our considerations, which pertaijH|-to-tbe| 
common experimental situation,riri Ep,Ep 3> max(e, IA5-1-I), we can employ the Andreev approximationLfcSSEfl 



and write kt 



kps 



k'p. It then follows that the appropriate wavevectors 



.5 ^ .,^^ = {2m/n^f/^[EF + PsKV'\ kf ^ .p. 
for the transmission of ELQ's and HLQ's are at angles d'g, ^S'g, with the interface normal, respectively. Within 
this approximation the components of the vectors k^ , k^ normal and parallel to the interface, can be expressed as 
kg = [ks, ^||s), and kg = {k'g, k^^g), in the F and S regions. From the conservation of fc[|5, we have then an analogue 
of Snell's law 



ips sm ( 



kpg sin 9-g, 



(2.2a) 



kps sin 6 = k'p sin O'g , 



(2.2b) 



which has several important implications, including the existence of critical anglesJlS as one encounters in well known 
phenomena in the propagation of electromagnetic wavesr^ 



Using the conservation of /ciig, the solution to Eq. (2.1), ^5 = {ug,v-^)'^, can be expressed in a separable form. 



effectively reducing the problem to a one-dimensional one. In the F region we write 



where 

ipsix) = 
analogously, in the S region we have 



Jksx 



*s(r) 



= ^»k||s-r 



i'six), 
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-ikgx 
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^s(.x), 



(2.3) 



(2.4) 



(2.5) 



tp'g{x) = Cge' 



dge^^'^'s'' 



e*'^~(e-rJs~/2e)i 
{e + VLg^/2e)^- 



ie + ng+/2e)2 
e-^*+(e-f]s+/2e)^ 

with ^g± = (e^ — |As±p)2, and the appropriate boundary conditionsoEa at the F/S interface are 

2mH 



V's(0)=V^(0), d.,^Jg{0) - d^^j'giO) 



-^siO)- 



(2.6) 



(2.7) 



We pause next to discuss some implications of Eq. (2^) for the various scattering processes. In typical realizations 
of fcrromagnct/HTSC structures, the appropriate FWM corresponds to Lq < 1.Cj Consider first Lq — 1, i.e. Ep — E 'p. 

I ^ u^. fQj. g^jj g — I incoming electron. Then, at any incident angle, Eq. (2.2) is 

> '" 



If X > it follows that kp^^ < k'p < kp-^ 



satisfied so that k\\ will be conserved. In this case \9\ > \9',\ > \9j\, and all the corresponding wave vectors are real. For 
an S —1 incident electron at angle 16*1 > | sm~^ {k'p /kpi)\, a solution of Eq. ( ^.2b| ) for a real 9'^ no longer exist, one has 
a complex 01 .c3 The scattering problem does not have a solution with propagating wavevectors in the S region: there 
is total reflection. The wavevectors for ELQ and HLQ have purely imaginary components along the x-axis, while their 
components parallel to the interface are real. This-cprresponds to a surface (evanescent) wave, propagating along the 
interface and exponentially damped away from it.t3 An analogous, but physically more interesting, situation occurs 
for AR in the particular case where \9\ is smaller than the angl e of to tal reflection and satisfies 16*1 > | sm~^ {kpi/kp^)' 



This regime corresponds to kuj > kp^. In this case it is Eq. ( ^.2a ) that has no solution for real angles. This means 



III " '^Fl- 

that Andreev reflection as a propagating wave is impossible. From the condition, which follows from the Andreev 
approximation, fc2.-|- k"^ = k"^., we see that the component Ahr along the x axis must be purely imaginary,E2l while k..-r 

is still real. With these considerations we then find 






-i{e 



FT ' 



1,2 U/2 



(2.8) 



where we have expressed k-r in terms of quantities which are always real and which pertain to the F region only. As 
with total reflection, there is propagation only along the interface and an exponential decay away from it. This case 
differs from that of total reflection in that, since the evanescence affects only the Andreev reflected component, there 
may still be transmission across the junction. 

The above considerations apply a fortiori in the presence of FWM. For example, if we now consider Lq < 1, we can 
see by inspection of Eq. (2.2), that there can also be total reflection for an S =| incident electron, when kpi > k'p. 
This condition would imply the absence of imaginary kj for any incident angle and any exchange energy. 



Returning now to the b asic equations, we see that by solving for ^s{x), '>P'si^) ^^ Eq. (2.4), ( |2.6| ) with the boundary 
conditions given by Eq. (2/7), we can obtain the ampUtudes as, bs, cs and ds, S =t,i- For each spin, there is a 
sum rule, related to the conservation of probability, for the squares of the absolute values of the amplitudes. We can 
thus, in a way similar to what was done in Ref. 36, express the various quantities in terms of the amplitudes 05 and 
bs only. These amplitudes are given by 



as 



AtsLsT+e~'^^+ 



Uss+U-. 



+'^ss~ 



Vss-V^s+^+^-e'^^'- 



-<t>s+) ■ 



(2.9) 



i(0S--0S+) 



Vss+U-ss- - Uss-V-ss+^+^-^' 

Uss+U^. - Vss^V-^.^r+T^eA<t's—^s+) 



(2.10) 



where we have introduced the notation r± = (e — ^s±)/|^S±|, Ls = Lq cos 6'g / cos 9 , describing FWM, ts = 
ks/kp^ = (l + psX)i/2, t- = kj/kpa^ ^ {1-psXy^^cose^/cose, for % real, (-i[(l + X) sin^ 61 - (1 - X)]i/Vcos6', 
for kj imaginary, see Eq. ( p.8D ). The other abbreviations are defined as: U-g^^ = t-g + ws±, Vss± = ts — ws±, 
ws± = Ls ± 2iZ, Z = Zo/cos0, where Zq = raH/fikp is the interfacial barrier parameter, as defined above. The 
limits Zq — > and Zq — > 00 correspond to the extreme cases of a metallic point contact and the tunnel junction limit, 
respectively. _ . 

Given the above amplitudes, the results for the dimcnsionless differential conductanceEEl can be written down in the 
standard way by computing, as a function of the excitation energy arising from the application of a bias voltage, the 
ratio of the induced flux densities across the junetifan.to the corresponding incident flux density. One straightforwardly 
generalizes the methods used in previous workoEjCj to include now the effects of unconventional superconductivity, 
FWM, and net spin polarization, to obtain, 



G ^ G| + G| 






\bsn 



(2.11) 



where we introduceLJjhe probability Pg of a n incident electron having spin S*, related to the exchange energy as 
P5 = (1 + psX)/2^ In deriving Eq. ( 2.11 ), care has to be taken to properly include the flux factors, which are, 
at X > 0, different for t he inc ident and the Andreev reflected particle. The ratio of wavevectors in the second term 
on the right side of Eq. ( 2.11 ) results from the incident electron and the AR hole belonging to different spin bands. 
The quantity k-^ in that term is real, the case of imaginary k-g can o»ly contribute to G-\ indirectly, by modifying 
|6||. It can be shownuS from the. conservation of probability currentEa that such a contribution vanishes for the 
subgap conductjaHce (e < |A|±|).c3 It is, furthermore, possible to express the subgap conductance in terms of the AR 
amplitude only.E3 At AT = we recover the results of Ref. |^. The suppression of the conductance due to ordinary 
reflection at X 7^ has the same form as for the unpolarized case since the magnitude of the normal component of 
the wavevectors before and after ordinary reflection remains the same. 

We fo cus i n this work (see results in the next Section ) on t he co nductance spectrum of the charge current as given 
by Eq. (2.11), but the amplitudes 05, 65, given by Eq. (2.9), (2.1C) can be used to calculate many other quantities of 
interest, such as current-voltage characteristics, the spin current, and the spin conductanccEiJ We consider also here 
angularly averaged quantities and notice that Eq. (2.11) implies that the conductance vanishes for \9\ greater than 
the angle of total reflection (we recall that this angle is spin dependent). We define the angularly averaged (AA) 
conductance, (Gs), as 



(G5 



de cos eGs{9)/ / de cos 61, 
ns Jus 



(2.12) 



where Vis is limited by the angle of total reflection or by experimental setup. This form correctly reduces to that 
used in the previously investigated spin unpolarized situation.Q One may choose a different weight function in per- 
forming such, angular averages, depending on the specific experimental geometry and the strengths of the interfacial 



However, all expressions for angularly averaged results, obtained from different averaging methods, 



nm 



scattermg 

would still have a factor of (1 -I- jp-jagp — l^sp) in the kernel of integration, and would merely require numerical 
integration of the amplitudes we have already given here. 



III. RESULTS 



A. Conventional pair potentials 



We present our results in terms of the dimensionless differential conductance, plotted as a function of the dinien- 
sionless energy E = e/Ao. We concentrate on the region E ^ 1 since for larger bias various extrinsic effects, such as 
heating, tend to dominate the behavior of the measured conductance.Ell While our findings, and the analytic results 
from Section O, are valid for any value of the interfacial scattering, we focus on smaller values oi Zq, Zq < 1, where 
the novel effects of ferromagnetism on Andreev reflection, and consequently on the conductance, are more pronounced 
than in the tunneling limit, Zq 3> 1. This regime on which we focus is also that which is believed to corresrinnd to 
several ongoing experiments of F/S structures, where the samples typically have small interface resist.aiji.ce.EjE3 To 
present numerical results, we choose E'p/Aq = 12.5, consistent with optimally doped FBa2Cu307_5.E3C3 We will 
include FWM, as parametrized by the quantity Lq introduced above, Ep = E'p/Lq. 

We first give some results for an s-wave PP, with a constant energy gap. This will serve to illustrate the influence 
of FWM coupled with that of Za within a simpler and more familiar context. In this case, for any incident angle, 9, 
of an injected electron the ELQ and HLQ feel the same PP with As± = Ag, and (j)s± = 0. Therefore, the results 
that we give here for the s- wave case and normal incidence (6 = 0), also correspond to the case of a PP of the dx^^y2 
form, with the angle a S (— 7r/2,7r/2), between the crystallographic a-axis and the interface normal, set to a = 0. 
This would represent an F/S interface along the (100) plane. 
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FIG. 1. G{E) (Eq.(2.11)) versus E = e/Ao. Results are for ^ = (normal incidence). The curves are for Zo = (no barrier): 
in panel (a) at exchange energy X = ho/Ep ~ (no spin polarization) they are (from top to bottom at any E) for the FWM 
values of Lq = Ep/Ep = 1, l/\/2, 1/2, 1/4, 1/9, 1/16. In panel (b) they are for X = 0.866. Since the curves now cross at E — 1 
they are drawn in different ways for clarity. For E > 1 they are in the same order as in panel (a) and for the same values of 
I/O, while for E < 1 they correspond, from top to bottom, to Lq ~ 1/2, l/\/2, 1, 1/9, 1/16. The Lq — 1/4 curve overlaps with 
that for Lq — 1 in this range. 



In Fig. ^ we show results for G{E), given by Eq. ( 2.11 ), at = 0, and Zq = (this limit of no interfacial barrier 
was also considered in Ref. pa). We plot results for various values of the FWM parameter Lq. Panel (a) corresponds 
to no polarization {X ~ 0) and panel (b) to high polarizati on X = \/3/2 w 0.866. For normal incidence, we have 
ts = (1 + psX)^/^, t-g~(\ — psX)^!"^ (as defined below Eq. (p.lC|)), and the subgap conductance can be expressed as 



G = 



32L^(1-X^) 



2U/2 



\hh + m + <i)io -^Ll- (t^t^ - m + ty)L^ + L2)r+r_|= 



(3.1) 



Panel (a) displays results in the absence of exchange energy. With increasing FWM (i.e. decreasing Lo)j the amplitude 
at zero bias voltage (AZB) decreases monotonically. This effect was explainedtZI in previous work as resulting from 
the increase in a single parameter .^e//, which combined Zq with the effects of FWM. Our curves with FWM (Lq < 1) 



reduce in the appropriate limits to those previously foundO with Lg = 1 and Zq -^ Zeff, Z^jf > Zq- We will see 
below that this is not the case at X 7^ 0. In panel (b) we give results for high X while keeping the other parameters at 
the same values as in panel (a). We notice that the presence of exchange energy gives rise to non-monotonic behavior 
in the AZB. At low bias, the conductance can be enhanced with increasing FWM (compare, for example, the Lq = 1 
and Lq = l/\/2 results), and form a zero bias conductance peak (ZBCP.) This behavior is qualitatively different 
from that found in the unpolarized case and the effect of |E,WM can no longer be reproduced by simply increasing 
the interface scattering parameter. Thus the often impliedll3'E2l expectation that the effects of Zq and Lq could also 
be subsumed in a single parameter in the spin polarized case is not fulfilled. In this panel we have an example of 
coinciding subgap conductances fo r Lq = 1 and Lq = 1/2. The condition for this coincidence to take place at fixed X 
can be simply obtained from Eq. (3.1) as 



M'^i/Lq - l'q 



{l-X^) = L'i{LQ^l) 



(3.2) 



where Lq, Lq correspond to two different values of FWM for which the subgap conductances will coincide. 

We next look, in the same situation as in the previous figure, at the effects of the presence of an interfacial barrier. 
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FIG. 2. G{E) for 6^ = and interfacial barrier strength Zq = 1. All the other parameters are taken as in the previous figure. 
In both panels, curves from top to bottom correspond to decreasing values of Lq. 



In Fig. g, we choose Zq = 1, while keeping all the other parameters the same as in the corresponding panel of 
the previous figure. In panel (a) we show results in the absence of spin polarization. A finite bias conductance peak 
(FBCP) appears at the gap edge. It becomes increasingly narrow with greater FWM (smaller Lq). Its amplitude is 2, 
independent of Lq. In panel (b), at X = 0.866 , the c onduc tance curves display similar behavior, but with a reduced 
FBCP at the gap edge. From Eqs. {1^, (|2.10|) , and (|2.1l| ), the amplitude of the FBCP in this case is 



G{E = 1) = 



4(^_ ^2)1/2 
l + (l-X2)l/2- 



(3.3) 



An interesting feature of this result is that it depends only on the exchange energy (spin polarization) and not on the 
FWM parameter or the barrier strength. It can be shown that this property holds for all angles of incidence. This 
is in contrast with the value of the zero bias conductance which depends, both for normal incidence and for other 
angles, on the value of the EWM. This dependence could introduce difficulties in the accurate determination of spin 
polarization from the AZB.E3E3 The gap edge value is less susceptible to these problems. 
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FIG. 3. Evolution of the zero bias conductance, G{E) for 6 — 0. Results are given at Zq — 1 for X determined from Eq. ( p.4| ) 
and values of Lq as in Fig. hi From top to bottom the curves correspond to values of {Lq,X) given by (1,0), (l/\/2, l/\/2), 
(1/2,0.866), (1/4,0.968), (1/9,0.994), and (1/16,0.998). 

The presence of spin polarized carriers, due to nonvanishing exchange energy, is usually heldE3c3^Ej to result in. 
the suppression of Andreev reflection and thus in a reduction of the subgap conductance. A simple cxplanation,E3 
which neglects the effects of FWM, predicts that the AZB should monotonically decrease with increasing X, because 
of the reduction of Andreev reflection, when only a fraction of injected electrons from the majority spin band can be 
reflected as holes belonging to the minority spin band. This follows from the reduction of the density of states in the 
minority spin band with increasing X, and eventually causes the subgap conductance to vanish for a half-metallic 
ferromagnet when X —* 1. 
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(a) 
X=0.866, Z 0=0 
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X=0.866, Zo=l 
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FIG. 4. {G{E)), the 6 averaged conductance, for an s-wave PP and the same values of X, Lo as in panels (b) of Figs, hi, Q, 
respectively. In both panels curves from top to bottom, ai E = 2, correspond to decreasing Lq. 

We now proceed to examine whether these findings are modified when FWM is taken into account. In Fig. y, 
which shows results at Zq = and normal incidence, we consider the evolution of the conductance curves for different 
values of X and Lq chosen to yield maximum AZB, {G{E = 0) = 2), starting from the step-like feature at Xp = 1 
and X = (see Fig. (m)). The condition for maximum AZB at fixed FWM and polarization can be derivedE3 from 
Eq. (O) and is 



fc|fc| = k'j: 



(i-x2)i/2=i; 



(3.4) 



We have used this equation to determine the optimal value of X for each value of Lq used in Figs. |l|, g. The 



resulting curves are plotted in Fig. ^ This figure reveals several interesting features. With the increa se of FWM and 
the correspondingly larger optimal spin polarization (according to the value of X found from Eq. (3^)), a ZBCP 
forms. This is a novel effect in which the peak arises from a mechanism completely different from the one usually put 
forward, where the ZBCP is attributed to the presence of unconventional superconductivity. In thatjCase, the ZBCP 
is produced by the sign change of the PP and the concomitant formation of Andreev bound states.ErQO Furthermore, 
if we compare these curves with those in panel (b) of Fig. |l|, we sec that the subgap conductance can increase with 
increasing spin polarization at fixed Lq- This implies that Andreev reflection can be enhanced by spin polarization. 

We now turn to angular averages (AA). In Fig. ^ we show angularly averaged results, obtained from the expression 
for (G), Eq. (2.12). The averaged results are no longer equivalent (as in the previous figures with normal incidence) 
to the case of a dx2_y2 PP with an F/S interface along the (100) plane: the angular dependence of the PP would 
then modify the results. Each of the two panels shown includes results for the same set of parameter values used in 
panels (b) of Figs. H, and 0, respectively. In panel (a) of the current figure we show how the novel features previously 
discussed are largely preserved after angular averaging. There is still formation of a ZBCP with increased FWM and 
the AZB retains its non-monotonic behavior with Lq, as in the case of fixed normal incidence. The angularly averaged 
results in panel (b), at Zq — 1, display behavior similar to that found in the 9 = case, with the conductance peak 
at E' = 1 becoming sharper at increasing FWM. 



B. Unconventional pair potentials 



We next consider an angularly dependent PP, specifically that for a d^2_y2 pairing state. With this PP we have 
different, spin dependent, PP's for ELQ's and HLQ's. These are given respectively by As± = Aq cos{29'g^), where 9'g^ 
can be expressed as 9'^^. = 9'g^a (we recal l tha t a is the angle between the interface normal and the crystallographic 
a- axis, and 9'g is related to 9 through Eq. (2.2)). 
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FIG. 5. G{E) for = tt/IO, a = tt/4, Zq = 0, and Lo = 1- In (a) the curves are for X = 0,0.5,0.7,0.8,0.866,0.95, (top to 
bottom at _B = 0). In (b), we plot the spin resolved conductance, for two values of X. The upper curve at _E > 1 corresponds 
to Gt, and and lower curve to Gi- 

In Fig. we give some of our results for d-wave pairing and a = tt/4 (interface in the (110) plane), in the absence 
of both interfacial barrier and FWM and at a fixed 9 = tt/IO, for various values of X. Panel (a) shows curves for 
the total conductance as it evolves from a step-like feature at X = to a zero bias conductance dip (ZBCD) for 
large spin polarization. The width of the plateau at X = is determined by a single energy scale set by the equal 
magnitudes of the PP's for ELQ and HLQ in that case, as given by As+ = As^ < A p, S =tji- As the exchange 
energy is increased, kp^ and fc^x are no longer equal. As one can see from Eq. (2.2b), it follows that 9i ^ 9\ and 



thus A|± ^ ^i±- These two different energy scales are responsible for the position of various features, such as the 
several finite bias conductance peaks (FBCP's) that are seen. 
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FIG. 6. G(E) for Q — tt/IO, a = 7r/6, Zo = 0, and Lo — 1. In both panels ordering and values of X for each curve are as in 
Fig. |. 

In panel (b) we show the spin decomposition G = G-^ + Gi, which better reveals these scales, at two different 
exchange energies. At X = 0.5, the shapes of G|, Gi are only shghtly modified from those in the unpolarized case. 
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FIG. 7. G{E) for 6 = vr/lO, Zq — 1, and Lo — 1. In both panels curves (top to bottom at £ = 0) correspond to 
X = 0, 0.5, 0.7, 0.8, 0.9. In (a) a = 7r/4, and in (b) a = 7r/6. 



At larger exchange energy, X — 0.866, the situation is very different, as shown in the figure. We also see, in panel 



(b), that as stated in the previous section, the evanescent wave associated with the imaginary kj does not contribute 
to the subgap conductance G| . 

In general, for an arbitrary orientation of the F/S interface, a ^ 0, 7r/4, at a fixed 9, aU the four spin dependent 
PP's for ELQ and HLQ will have different magnitudes. There are, therefore, specific features at four different energy 
scales. It is only for the particular and atypical (but often chosen in theoretical work) case of a = 7r/4 that these four 
scales reduce to two. 

In Fig. owe show the general behavior by choosing a — 7r/6, while retaining the values of all the other parameters 
from the previous figure. One can easily calculate, for example, that at X ~ 0.5 the normalized values of the PP are, 
in units of the gap maximum, Ao, |A|+| = 0.963, |A|+| = 0.250, |A|+| = 0.822, |Ax_| = 0.083. These numbers can 
also be approximately inferred from the spin resolved results given by the solid lines in panel (b). 
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FIG. 8. {G{E)}, at Zq = 1, and Lo = 1 for d^2_ 
As± = Aocos(2Sg^) + iO.lAo. In panel (a) a = 7r/4 and in (b) a = 
correspond to X = 0, 0.5, 0.9, in both panels and for each pair potential 



and d^2_y2 + is pair potentials. The latter is of the form 
7r/6. From top to bottom (at E = 0), the curves 



We next turn to the case where there is a nonvanishing potential barrier, choosing for illustration the value Zq — 
1. In the ^Jfgwce of spin polarization, the formation of a ZBCP at finite barrier strength has been extensively 
investigatedafl'Lj and explained by Andreev bound states in the context of d-wave superconductivity. We will consider 
here also the effects oi X , not included in previous work. In Fig. ^ we show results for various values oi X at a — 7r/4. 
(in panel (a)) and a — 7r/6 (panel (b)). One can see that for intermediate values of X the conductance maximum is 
at finite bias. Comparing the two panels, one sees that the AZB at a fixed X ^ is larger for a — 7r/4, in agreement 
with the results obtained for the unpolarized case where, at zero bias, the spectral weight is maximaB for a (110) 
interface. For a different choice of incident angle 9 there will be, if the values of all other parameters are held fixed, 
a change in the effective barrier strength for various scattering processes. We recall (see below, Eq. ( |2.10D ), that 
Z = Zq/cos9, and with an increase in 16*1 typically there will be, as in the unpolarized casej^ a decrease in the 
amplitude for Andreev refiection and an increased amplitude for ordinary refiection. 

Results such as tho se di scussed above can be obtained as a function of angle, and the angular average can then be 
computed from Eq. ( 2.12 ). We will combine showing some of these angularly averaged results with a brief study of 
another point: it is straightforward to use the formalism discussed here to examine more complicated superconducting 
order parameters. A question that has given rise to a considerable amount of discussion is that of whether the 
superconducting order parameter in high Tc materials is pure d-wave or contains a mixture of s wave as well, with 
an imaginary component, so that there would not be, strictly speaking, gap nodes, but only very deep minima. With 
this in mind, the effect of a possible jlUmaginary" PP admixture (for example in a d + is form) on Andreev bound 
states has also been recently studied.lljilj We consider this question here, including the effects of polarization. In 
Fig. H, we illustrate the difference in the angularly averaged conductance values obtained for a pure (1^2 _y2 PP and 
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for a mixed d^'^-y'^ + is case. We choose the particular form Ag^ = 0.9Ao cos(20^_|_) + iO.lAo. The phase of the PP, 
4's±: is no longer equal to tt or as in the pure d-wave case. We give AA results for several values of X, both for the 
pure d and the mixed d-^is cases. 
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FIG. 9. G{E) for e = n/10, Zq = 0, and Lo = 1/2, Lq = 1/3. In panel (a) a = 7r/4 and in (b) a = 7r/6. From top to bottom, 
at _B = 0, curves correspond to X = 0, 0.5, 0.8, in both panels and for each pairing potential. 

The former represents the angular average of results similar to those previously displayed. 
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FIG. 10. Conductance curves for 9 — tt/IO aX Zq — \ with the same parameters and ordering as in Fig. H. 

As in the unpolarized case,oEj the is admixture in the PP is responsible for a FBCP, approximately at E — 0.1. 
The conductance maximum is reduced with increased X and with departure from a (110) oriented interface. Replacing 
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the dj^2_y2 + is PP by a "real" admixture d^2_y2 + s (taking again O.IAq for the s-wave part) gives results almost 
indistinguishable from the pure d-wave for any value of spin polarization. 

To show the effects of FWM on conductance for a pure d-wave PP we take Lq = 1/2, 1/3 and give results at the 
fixed angle, 9 = tt/IO, previously considered. In Fig. ^ wc show curves at Zq — and a — n/4 (panel (a)), and for 
a = 7r/6 in panel (b). It is useful to compare this figure to panel (a) in Figs. H, and 0, corresponding to no FWM 
for a — 7r/4 and 7r/6, respectively. In the absence of spin polarization the effect of FWM resembles the influence of 
a nonvanishing barrier strength, Zq, and leads to the formation of a ZBCP, which becomes increasingly narrow for 
smaller Lq. The effect of moderate spin polarization (X < 0.5, for comparison with the above mentioned figures) on 
the AZB is rather small for Lq = 1, 1/2 but it is significantly larger at Lq = 1/3. In the next figure. Fig. nO, we use 
Zq = 1 and the same parameters as in the previous figure, so that the infiuence of barrier strength can be gauged. 
One sees that in the presence of spin polarization the position of the conductance maximum depends on FWM. With 
increasing mismatch, the FBCP evolves into a ZBCP. By comparing the curves corresponding to Lq = 1 in Fig- Q 
with those for smaller Lq in Fig. 11^ , it is interesting to notice that an effect similar to that discussed previously 
for s-wave PP without an interfacial barrier and at normal incidence is also manifested in other regimes, in that the 
conductance maximum can actually be enhanced, in the spin polarized case (at fixed X), by the FWM. 

IV. CONCLUSIONS 

In this paper we have studied the conductance spectra of ferromagnetic/superconductor structures. The expressions 
for Andreev reflection and ordinary reflection amplitudes which we have given, allow one to simply ahtain other 
quantities of interest such as current- voltage characteristics or conductance spectra for spin current EJ We have 
developed the appropriate extensions of the standard approach and approximations used in the absence of spin 
polarization. This has enabled us to present analytic results. Within these approximations, and with the inclusion 
of FWM, we have shown a number of important qualitative differences from the unpolarized case or from that where 
spin polarization is included in the absence of FWM. r-ir-l 

Our considerations may also be important in the interpretation of recent experimentsli3i23 attempting to use tun- 
neling to measure the degree of spin polarization in the ferromagnetic side of the junction, since the experimental 
determination of spin polarization in a ferromagnet is a very difficult and important experimental question in its own 
right. As we have shown, the ZBCP is sensitive to both spin polarization and FWM, while the gap edge amplitude 
depends only on X. It is then not possible to straightforwardly determine the spin polarization by using the results for 
the amplitude of the zero bias conductance unless the appropriate FWM of the E7S structure is known and properly 
taken into account. Furthermore, FWM can not, unlike in the unpolarized case,EZI be simply described by a rescaled 
value of the interfacial barrier strength. 

The procedures used here have the advantages of simplicity and of allowing for analytic solutions. These advantages 
have enabled us to investigate widely the relevant parameter space. We have left for future work considerations that 
would have diminished these advantages. Among these are the question of the self consistent treatment of the PP, 
inclusion of spin-fiip scattering or of a more realistic band structure, and non-equilibrium transport. However, we 
believe that the the methods we have employed are sufficient to elucidate the hitherto unappreciated subtleties and 
the richness and variety of the phenomena associated with spin polarized tunneling spectroscopy. 

We hope that our work, as reported here and in Ref. ^, will prompt additional experiments and theoretical work. In 
particular, an important clue about spin polarized transport would be provided by measurements of the spin resolved 
conjductance. Indeed, we have already become aware of two very recent related preprints leading intOi-these directions, 
onco on Andreev reflection and spin injection into s- and d-wave superconductors, and anotherES discussing, in 
conventional superconductors, out of equilibrium enhanced Andreev reflection with spin polarization. 
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